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A new transient effective theory of the relativistic Boltzmann equation is derived for locally 
momentum-anisotropic systems. In the expansion of the distribution function around a local “quasi¬ 
equilibrium” state a non-hydrodynamic dynamical degree of freedom is introduced at leading order 
that breaks local momentum isotropy. By replacing the deviation of the distribution function from 
this quasi-equilibrium state in terms of moments of the leading-order distribution and applying a 
systematic power counting scheme that orders the non-hydrodynamic modes by their microscopic 
time scales, a closed set of equations for the dynamical degrees of freedom is obtained. Truncating 
this set at the level of the slowest non-hydroynamic mode we find that it exhibits transient oscillatory 
behavior - a phenomenon previously found only in strongly coupled theories, where it appears to be 
generic. In weakly coupled systems described by the Boltzmann equation, these transient oscillations 
depend on the breaking of local momentum isotropy being treated non-perturbatively at leading 
order in the expansion of the distribution function. 

PACS numbers: 


I. INTRODUCTION 

The Boltzmann equation provides a microscopic link 
between statistical physics and macroscopic fluid dynam¬ 
ics. For weakly coupled, dilute plasmas of massive parti¬ 
cles for which two-particle correlations can be neglected 
the Boltzmann equation with only elastic binary colli¬ 
sions yields a good description of the microscopic physics. 
While hydrodynamics provides an effective theory for the 
evolution of conserved macroscopic quantities (energy, 
momentum, and conserved charges), a full description 
of the microscopic dynamics of the Boltzmann equation 
must include an infinite set of non-hydrodynamic modes. 
Hydrodynamic modes have dispersion relations that sat¬ 
isfy linifc^o w„(fc) = 0; all other modes are called non¬ 
hydrodynamic. The former dominate the dynamics of 
the system at large spatial and temporal length scales. 
The influence of non-hydrodynamic modes is suppressed 
for small values of the Knudsen number Kn (the ratio 
between the mean free path and the typical macroscopic 
length scale of the system). 

Systematic coarse-graining leads to effective theories 
with varying degrees of freedom and non-hydrodynamic 
modes. For relativistic systems, causal effective theories 
must include at least the first non-hydrodynamic mode 
in the dynamics of the dissipative currents [1, 2]. Start¬ 
ing from the Boltzmann equation and expanding the dis¬ 
tribution function around an isotropic local equilibrium 
state it was shown [1] that all non-hydrodynamic modes 
have frequencies that for k —> 0 lie on the negative imagi¬ 
nary axis of the complex to plane. At large times, only the 
slowest non-hydrodynamic mode with the smallest fre¬ 
quency |w| survives. This, together with a power count¬ 
ing scheme in Knudsen number Kn and inverse Reynolds 
numbers R” 1 (the ratios between the dissipative currents 
and the thermal equilibrium quantities), allowed Deni- 
col etal. to systematically derive a resummed transient 
relativistic fluid dynamical theory [3]. 


The approach in [3] is constrained by expanding the 
Boltzmann equation around an isotropic local equilib¬ 
rium distribution, and hence by treating all deviations 
from local momentum anisotropy perturbatively. The 
infinite moment hierarchy is truncated by a power count¬ 
ing scheme based on inverse Reynolds numbers that in¬ 
clude all dynamically generated momentum anisotropies. 
When a system expands strongly anisotropically, such as 
the quark-gluon plasma fireball created during the early 
stages of a relativistic heavy-ion collision [4, 5] , these lo¬ 
cal momentum anisotropies can become large and must 
be treated nonperturbatively. Here we do so by modify¬ 
ing the local thermal equilibrium distribution by intro¬ 
ducing a non-hydrodynamic degree of freedom describing 
the deviation from momentum isotropy, thereby defining 
an anisotropic quasi-equilibrium state [6] . In a procedure 
that closely follows [3] , but now resums and truncates the 
moment hierarchy according to a modified power count¬ 
ing scheme based on the Knudsen and residual Reynolds 
numbers from which the largest contributions arising 
from the local momentum anisotropy have been elimi¬ 
nated [7], we derive a new type of transient relativistic 
fluid dynamics in which the slowest non-hydrodynamic 
mode (associated with the local momentum anisotropy) 
turns out to exhibit transient oscillations. 

Such oscillations appear to be generic in macroscopic 
effective theories of strongly coupled plasmas [2, 8-11], 
but were thought not to exist in analogous theories of 
weakly coupled systems that are based on the Boltzmann 
equation [1, 2]. We find that they arise generically also 
in weakly coupled systems and are absent only if the 
macroscopic effective theory is based on a perturbative 
expansion around an isotropic local equilibrium state. 
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II. DERIVATION 

We start from the relativistic Boltzmann equation for 
a gas of particles with mass to, 

k^d,jf(x,k) = C[f](x,k), (1) 

where k^ = (Ek,k), with E k =\ / k' 2 +m 2 , is the particle 
four-momentum. We write shortly /*. = /(:£, fc), Ck = 
C(x,k), and use dK = d 3 k/[(2TT 3 )Ek] to denote the 
Lorentz invariant momentum integration measure. The 
collision term for elastic 2-body collisions is [12] 

C k [f\ = \ f dK'dPdP'W kk '^ PP ’ 

x (/p/p'/fc/fc'-/fc/fc'/ p /p') , ( 2 ) 

where Wkk'^pp' is the transition rate for the process k + 
k' —> pAp' and / = 1 — af 1 where a = ±1,0 for fermions, 
bosons, or classical distinguishable particles, respectively. 

We consider systems undergoing highly anisotropic col¬ 
lective expansion, assumed to be stronger in the lon¬ 
gitudinal (z) than in the transverse (r± = (x,y)) direc¬ 
tions. In dissipative systems this leads to a local mo¬ 
mentum anisotropy which, for simplicity, we take at 
leading order to be azimuthally symmetric around z in 
the local rest frame (LRF), treating non-leading cor¬ 
rections to the momentum anisotropy as small pertur¬ 
bations. These assumptions are implemented by writ¬ 
ing [13] f k = f a ,k+Sfk = /a,fc(l+/a,fe0fe) where in 
the LRF the momentum-anisotropic (“a”) leading-order 
term takes the form [14] 

/a.fc = fe q {Pa.\Jm 2 +k\+{ l+£)^ ) ( 3 ) 

with the equilibrium distribution / eq (C) = l/(e^+a). The 
deviation <f>k is assumed to be <1. £(x) parametrizes 
the leading-order momentum anisotropy at point x, and 
/3 a (x) is the local “transverse temperature” at that point. 

The dynamics of the Boltzmann equation is character¬ 
ized by a hierarchy of microscopic time scales. The tem¬ 
perature parameter /3 a , related to the energy density at 
point x whose evolution is controlled by energy conserva¬ 
tion, varies on hydrodynamic time scales. The evolution 
of the anisotropy parameter £, like that of all other dissi¬ 
pative flows, is not constrained by conservation laws and 
thus happens on the microscopic time scales associated 
with the Boltzmann collision term. We are interested in 
the “transient anisotropic regime” where the fast micro¬ 
scopic processes have all decayed except for the ones asso¬ 
ciated with the slowest of these microscopic time scales. 

The conserved particle current and energy-momentum 
tensor are given by the first and second moments of the 
distribution function (with (■••) = f dK (• • •) /*,): 

= (jfc' 1 ), = {Wk v ). (4) 

The velocity held u^ (Landau frame) is defined as the 
timelike eigenvector of T^: T^u v = Eu U Its eigenvalue 


£ is the LRF energy density. We decompose the four- 
momenta into temporal and spatial LRF components: 
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W = ( u-k K + k W = E k u* + K i X i■ (5) 

i=1 

The four-vectors X? reduce in the LRF to unit vectors 
along the 3 spatial directions, and (Ek,K) are the LRF 
components of Using the generalized Landau match¬ 
ing condition that /3 a in (3) be chosen such that 5f k does 
not contribute to the LRF energy density £, one finds the 
anisotropic viscous hydrodynamic decomposition [13] 

= Af u » + ft* , 

r = W-(?x + + P'**'+ 7T # “', (6) 

Af=(u-k) a , Af^ = (k^)~ s , 

S = ((u-k) 2 ) a , fl = -i(A a Pk a k p ) s , m 

7T^=<fc<^>)j, V L =(k 2 z ) a , U 

V± = \((kl + k 2 y )) a , = (V L - V±) z *z v , 

with (•••) a = / dK(---)f a>k , (•■■)!= J dK(- ■ ■ )Sf k , z» 
= Ag (i.e. the z unit vector in the LRF), and A^ v '> = 
Awhere A^ = A&* A” ] - A^A aP /3 and A^ = 
g^-u^uE = ~T,i X i X i ■ Note that A M defined in (5) 
can also be written as A^ = A^A u . 

Following [3], we now replace the Boltzmann equation 
by a hierarchy of moment equations. We expand the devi¬ 
ation <j>k (x) in a complete orthogonal basis of irreducible 
tensors 1, k^\ k^k ^, k^k v k^\ ..., constructed from 
the spatial components of k in the LRF at point x [15]: 

oo Ni 

Mx) =^(^c^ 1 --^ e) (x)P^ ) (x,k- u (x))^k {tll ---k tle) . 

£=0 n =0 

( 8 ) 

The Pn ^ = Ylr=o a nr (x)[u(x)-k) r are orthogonal poly¬ 
nomials of order n in the LRF energy u-k. Their orthog¬ 
onality can be used to express the expansion coefficients 
Cn 1 " ^ (x) in terms of the irreducible moments of the 
residual distribution Sfk [3, 16] 

p(W"^( x) = J dK(u-k) r k ^ ■ ■ ■ k^Sfk ■ (9) 

After some algebra [17] the distribution function is re¬ 
written as a series in the irreducible moments: 

oo oo iVfc 

h = /a, fc (i+ f\ k EEE (< fc) ) $::!% {x, u-k) 

l —0 k =0 r=0 

x ■■■k w) ) , (10) 

-her e 
with 

«>)£:$(*)= CD 

m r. 

E / dK(u-k) r+s k ^ ■ ■ ■ k^k (Vl ■ ■ ■ K k) f a ,kf a ,k- 

s—0 J 
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From here on we consider for simplicity a gas of mass¬ 
less particles with vanishing chemical potential [17] for 
which energy and momentum are the only conserved 
quantities. Using (6), (7), and (3), their conservation 
law d^T^ 1 ' = 0 then yields the following hydrodynamic 
equations (see [13] for the definition of the modified ther¬ 
modynamic ^-integrals over f a . k ): 


/3a = ‘ ^ 


J2,0, —1 


(£ + 


2V ±+Vl \ 

3 J 
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-(v^+a -Ht 


fijy 


( 12 ) 


u" = - A^ Q (P Q/3 +7f a ^)] I ( £+V ±) . 


While the scalar and vector irreducible moments, p r 
and p(f, control the bulk viscous pressure and heat cur¬ 
rent (which vanish here), we assume that they will not 
play important roles in the shear dynamics. We will also 
ignore contributions from irreducible moments of rank 
> 2 since they are of higher order in the power-counting 
scheme discussed below [3]. The dynamical equations 
of motion for p£ v are obtained [3, 18] by applying the 
operator where D = u-d is the time derivative 

in the LRF, to its kinetic definition (9) and using the 
Boltzmann equation (1) in the form (/ = Df ) 

Sfk = - A,fc - X (k-V(U, h +6f k ) - c fc [/]) . (13) 


Here is the velocity shear tensor and 6 = Linearizing the collision term in the deviation <p k one 

d-u the scalar expansion rate. finds [17, 18] 


n 2 


KVpf + Y.^rXyPn^ = W V i + M? X Z A + (affro + {°®Y VXP * Ap + 


n =0 


z I 

+^r J —^ ( r a p - -(2r+5)p^a? + 2 p^X -(r+4 )% v 0. (14) 

J 2,0,-1 ' 6 


( 2^ (2\ (2~\ 

Here w pl/ = (X p u v — V„Up)/2 is the vorticity tensor. The tensors C, A4, B, a g , a), , and aL are built from zi M , 
z^ and A**", multiplied by thermodynamic integrals over the anisotropic leading-order distribution f ak that depend 
on £ and /3 a ; their explicit expressions will be given in [17]. Our notation follows the convention used in [3] as closely 
as possible. In (14) all comoving time derivatives of /3 a and u^ were eliminated in terms of spatial gradients and 
comoving time derivatives of £ by using Eqs. (12). The matrix A^ 2 \ which couples moments of different orders r, 
embodies all the microscopic information contained in the linearized Boltzmann collision term: 



1 

2 


dKdK'dP dP' W kk ^ pp , 


X 



(n 


(2j2) \ (Act) 7 / ,, 
k'n r<r) 


(15) 


The main step to close this system of equations is to 
identify and separate the microscopic time scales. This is 
accomplished by determining the eigenmodes of A ^. Its 
dimension, which determines the number of eigenmodes, 
is given by the cutoff N/_ in the order of the polynomi¬ 
als Pn' 1 in the expansion (8) of the deviation (j) k which 
controls the accuracy with which the dependence of cf> k 
on the LRF energy can be resolved. For N( —> 00 one 
obtains an infinite number of eigenmodes and associated 
microscopic time scales. 

Introducing the matrix fthat diagonalizes A^ 2 \ 



multiplying (14) by ((f l and summing over r 

yields 


xX + ixfXpXf = x+{a$r v 0 (16) 

+{a^y vXP(J *p + ( Xi) IJ,, ' XPu} \p + higher-order terms, 


where Xp v = £^! 0 ((n XjXp X are the eigenmodes, 

and ^r Ai - A " = eSo (X^)7xf Xl ''' Xn for ^ e 

{C , M , <Yg\ aY\ a^}- “Higher-order terms” arise from 
those proportional to p£ v on the r.h.s. of Eq. (14). 

Equation (16) shows that the elements of the inverses 
of the “eigenvalue tensors” [x^Ya'p can be interpreted 
as microscopic relaxation times. Multiplying (16) by the 
inverse tensor [(x -1 )?/]^ and taking it to zero while 
keeping [{x~ l X]YYp etc. fixed, we see that the term 

involving x'f"^ is also driven to zero and that the eigen¬ 
mode assumes its “Navier-Stokes” limit: 

xr * + Kx-XYcYpM^z A (17) 

+[UXYX(Xr ] ) a0 o + [(xX?Tp(XX 0Xp ^ P 

+[(x _1 )[ 2) ]^(dLr) a/3Ap WAp + higher-order terms. 
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The terms explicitly listed on the r.h.s. of Eq. (17) are 
proportional to the Knudsen number Kn; the unlisted 
“higher-order terms” all involve an additional factor pro¬ 
portional to one of the modified inverse Reynolds num¬ 
bers ~ [19]. 

Next we order the elements of the eigenvalue tensors 
[Xr 2) ]^ by magnitude such that [xr 2) ]^ < [Xr+i]^- 
We consider the transient regime when the slowest mode 
(r = 0) still evolves dynamically according to Eq. (16) 
while all faster modes with r > 1 have already reached 
their asymptotic Navier-Stokes limits (17). Inverting the 
relation between eigenmodes and irreducible moments, 
(>T = Yhf=o(^iT)apXjP, normalizing the diagonalizing 
matrix by requiring (^oo^)a /3 = 1-Aand remembering 
f -T^=p^ v , Eq. (17) yields 

pr = (*4 2) )^+ fft v i +^r A iA + dfre 

+ {€iY VXp °Xp + {€iY vXp “x P + h- 0 .terms. (18) 

Here we defined i in terms of C by 

n 2 

r^E (^ 2) )%£?, (19) 

r—0 


and similarly for m, T)g 2 \ fffl and f in terms of A4, 
oi'Pj offl and aS\ with the microscopic relaxation times 

n 2 

= E^^kx- 1 )^]^^- 1 )^]^;. ( 20 ) 

m—0 


Equation (18), built on Eq. (17), implements the Navier- 
Stokes limit for all microscopic eigenmodes except the 
slowest one, r = 0. It expresses all higher-order tensor 
moments in terms of the lowest-order one, Pq U = , 

which appears in the energy-momentum conservation 
law. This is the key step in closing the equation of mo¬ 
tion (14) to simplify the coupling to higher order irre¬ 
ducible moments. To this end we multiply Eq. (14) by 
(tw)^ = [(^l (2) )“ r 1 ]^ and insert Eq. (18). Some alge¬ 
bra yields a set of equations, labeled by the order n of the 
irreducible moment considered, of the following form: 


+{v { f>y v + {v ( y ) y v 

= {T 1 { o' n) y v o + {r) > £' ri) y vat 3 
+(j ( ' n '>y i/ _|_ ( jc^y v . 


+ {r^y v i + (a ^y vX z x 
<T a p + {rj£' n) y vap w afi 

( 21 ) 


The coefficients of the various terms are related to 
(as indicated by the chosen notation) but not identi¬ 
cal with those in Eq. (18); their explicit expressions 
will be given in [17]. The first two terms on the left 
and the terms on the right of Eq. (21) are similar in 
structure to those of the generalized Israel-Stewart the¬ 
ory (DNMR theory) derived in [3]; the more complex 


Lorentz tensor structure results from the local momen¬ 
tum anisotropy of our leading-order distribution func¬ 
tion (3). As in [3], the J and K. terms are 0(Kn 2 ) and 
0(Kn • Rf 1 ), respectively [20]. The next two terms on 
the left of (21) arise from the second and third term 
on the right of Eq. (18), originating in the momen¬ 
tum anisotropy of the leading-order distribution / a *.. 
The last two terms on the left of Eq, (21) originate 
from time derivatives of these terms when evaluating 
the first term in Eq. (14). They include (in ('D^y 1 ') 
second time derivatives ~ £, z x of the parameters 
£, z x characterizing the deviation from local momentum 
isotropy in / a fc, as well as products of first time deriva¬ 
tives of these parameters with other first derivatives 
(~ z-z, z x zP , £i\ £-0(Kn, R” 1 ), i A -0(Kn,R t r 1 )) (these 

are collected in (D^y u ). Such terms occur gener- 
ically whenever the Boltzmann equation is expanded 
around a leading-order distribution function that is lo¬ 
cally anisotropic in momentum space [21]. As we will 
see below, the terms involving second time derivatives of 
£ identify this anisotropy parameter as a quasi-normal 
mode of the non-equilibrium dynamics which undergoes 
transient oscillations that are damped on a microscopic 
time scale that is related to the elements of (r ■^ > y v and a 

similar coefficient multiplying the £ term in (T>^y u [22], 
Such quasi-normal modes appear to be a generic feature 
in strongly-coupled theories which cannot be decribed by 
the Boltzmann equation [2, 8-11]. 

For n = 0, Eq. (21) controls the dynamical evolution 
of the residual shear stress fIt is coupled to the 
dynamical evolution of £ for which we can take Eq. (21) 
for n = 1. (When expanding around an isotropic leading- 
order distribution this equation is not needed.) Eqs. (21) 
for n= 1, together with its siblings for the bulk viscous 
pressure and heat flow vector, are sufficient to solve for 
all possible momentum-space deformation parameters of 
the leading-order distribution (9 at most). For n> 1, 
Eqs. (21) are only needed if we want to include dynamics 
at even shorter time scales, by making p p ^ dynamical 
variables instead of using their Navier-Stokes limits (17) 
as we have done here. 


III. EIGENMODE ANALYSIS 

To obtain dispersion relations for the eigenmodes of 
the linearized equations of motion, we consider a linear 
perturbation around a static background [15]: 

£ = £ 0 + e iut ~ ikz j£ ; 

Tf'"' = + e iut ~ ikz 8 ^ v , 

«** = < + e iult ~ ikz SuP , ( 22 ) 

with £o= const. 7 Tq I '= 0, and = (1,0,0,0), re¬ 
spectively. To study the shear channel we consider flow 
in the ^/-direction with a gradient in the z-direction, 
Su 11 = (0,0, Su v (t, z), 0), and study jr yz . For simplicity, 
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we assume 5z M = 0, no vorticity, and constant tempera¬ 
ture; the system is then not expanding, i.e. 9 = 0. In 
the linear regime, the elements of , (J7^")) A ‘ ! ', 

( K ,are all terms of 0(S 2 ) and can therefore be ig¬ 
nored [25]. In only one term ~ £ contributes 

[17] which we here denote as (^l"))/^. Introducing the 
shorthand notations 


ics at early and hydrodynamic behavior at later stages 
of a nuclear collision is characterized by transient os¬ 
cillations of the slowest non-hydrodynamic modes, with 
Im u) 


= (rji 2 ’ n) ) yzyz , 

N 

J 

1 

III 

J 

* -Uy> 




% 

III 

III 



' too 

O = (ri( 2 )\yz 
“n - V i6 n0 Jyz * 

d = d(P vz , 

(23) 

• 

• 


the linearized equations of motion can then be written as 
Ax = 0, with 


/iuj(£ + V '_ l ) — ikd —ik \ 

A = I — ikr — uj 2 £(°) + iw^ 0) iujri 0 ' 1 + f2 0 I , 

y — ikr /W —oj 2 ^ 1 ) + iuiri 1 ' 1 + VL\J 

x = (8u y , <$£, Stt vz ) T . (24) 

We note that the coefficients ^ n ' > vanish when the ex¬ 
pansion (10) is truncated at A ^2 = 0, so in Eq. (14) ATj 
must be > 1. For nontrivial solutions of (24) the deter¬ 
minant of A must vanish. In the long wavelength limit 
k —> 0 this determinant reduces to a fourth-order poly¬ 
nomial in cj( 0) = limfc_>.o w(A;), with two degenerate zeros 
at Wo = 0 describing hydrodynamic modes, and two non¬ 
hydrodynamic modes corresponding to solutions of the 
quadratic equation 

+ i w (r(»>r< 1 '-r('>T< 0 >+n 0 |(«-n l{ ») 

+ (si„T 5 (11 -n 1 T s ( ">) = 0 (25) 

that have the form u>± = ±u> r —icvi, with a; r , u>i > 0 (see 
Fig. 1). [That cjj > 0 is not obvious from (25) but neces¬ 
sary for the system to dynamically evolve towards local 
momentum isotropy.] The equal and opposite real and 
the negative imaginary parts of the frequencies of the 
non-hydrodynamic modes indicate that the momentum 
anisotropy parameter £ undergoes transient (damped) os¬ 
cillations with frequency uj r and damping rate uji. Similar 
transient oscillations were found in strongly coupled plas¬ 
mas when described using the fluid/gravity correspon¬ 
dence [2, 8-11]. 

IV. CONCLUSIONS 

In strongly coupled systems it was previously shown 
[11] that the transition between pre-equilibrium dynam- 


FIG. 1: Shear channel eigenmodes in the zero wavenumber 
limit. Shown are the two degenerate hydrodynamic modes at 
(Jo = 0 and the first two non-hydrodynamic modes at cj± that 
exhibit transient oscillatory behavior. 


(overdamped) oscillation amplitudes decaying on micro¬ 
scopic time scales. The present work, which is based on 
weakly coupled microscopic dynamics using the Boltz¬ 
mann equation, suggests that these oscillations arise 
from large deviations from local momentum isotropy in 
rapidly and anisotropically expanding systems: When 
strongly anisotropic expansion drives the local momen¬ 
tum distribution away from isotropy, the parameters 
describing this anisotropy do not decay exponentially, 
but perform damped oscillations. The matrix A in 
Eq. (24) couples their oscillations to Tr^ which there¬ 
fore oscillates in sync. These transient oscillations are a 
generic phenomenon, but capturing them in the macro¬ 
scopic equations of motion requires the inclusion of 
expansion-driven local momentum anisotropies in the mi¬ 
croscopic distribution function already at leading order. 
If the distribution function is instead expanded around 
a locally isotropic equilibrium state when deriving the 
macroscopic equations of motion from moments of the 
Boltzmann equation, only exponentially damped non¬ 
hydrodynamic modes are found [1]. 
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